SDSS J0903+5028: A New Gravitational Lens 



David E. Johnston, 1 ' 2 Gordon T. Richards, 3 Joshua A. Prieman, 1,2 ' 4 Charles R. Keeton, 1 ' 5 
Michael A. Strauss, 3 Robert H. Becker, 6 ' 7 Richard L. White , 8 Eric T. Johnson, 2 Zhaoming Ma, 1 
Mark SubbaRao, 1 ' 9 Neta A. Bahcall , 3 Mariangela Bernardi , 18 Jon Brinkmann, 10 Daniel J. 
Eisenstein, 11 Masataka Fukugita, 12 Patrick B. Hall, 3 ' 13 Naohisa Inada, Gillian R. Knapp, 3 
Bartosz Pindor, 3 David J. Schlegel 3 Ryan Scranton, 15 Erin S. Sheldon, 1 ' 2 Donald P. Schneider , 
C*~) ■ 16 Alexander S. Szalay, 17 and Donald G. York 1 

o 

ABSTRACT 

We report the discovery of a new gravitationally lensed quasar from the Sloan Dig- 
ital Sky Survey, SDSS J090334.92+502819.2. This object was targeted for SDSS spec- 
troscopy as a Luminous Red Galaxy (LRG), but manual examination of the spectrum 
showed the presence of a quasar at z ~ 3.6 in addition to a red galaxy at z = 0.388, and 
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the SDSS image showed a second possible QSO image nearby. Follow-up imaging and 
spectroscopy confirmed the lensing hypothesis. In images taken at the ARC 3.5-meter 
telescope, two quasars are separated by 2'.'8; the lensing galaxy is clearly seen and is 
blended with one of the quasar images. Spectroscopy taken at the Keck II telescope 
shows that the quasars have identical redshifts of z ~ 3.6 and both show the presence 
of the same broad absorption line-like troughs. We present simple lens models which 
account for the geometry and magnifications. The lens galaxy lies near two groups of 
galaxies and may be a part of them. The models suggest that the groups may contribute 
considerable shear and may have a strong effect on the lens configuration. 

Subject headings: gravitational lensing — quasars: individual (SDSS J090334. 92+502819. 2) 

1. Introduction 

Gravitational lenses have become important astrophysical and cosmological tools in several 
ways. The frequency of lensing is in principle sensitive to the dark energy density (Fukugita et al. 
1990, Turner 1990, Fukugita & Turner 1991, Kochanek 1995, but see Keeton 2002) and the mat- 
ter density of the Universe (Mortlock & Webster 2000a). Lens statistics also probe the properties 
of the lensing galaxy systems, such as their mass distribution, potential well depth and extinc- 
tion (Chen et al. 1995; Chae et al. 1998; Keeton et al. 1998; Malhotra et al. 1997; Keeton 2001a; 
Keeton & Madau 2001). In addition, measurement of the time delay between images in individual 
lensed quasars can be used to measure the Hubble parameter (Refsdal 1964). 

Since the discovery of the first double quasar Q0957+561 (Walsh et al. 1979; Schild & Thompson 1997), 
about 80 lensed quasars have been discovered 19 . In this paper, we report the discovery of another 
lensed quasar, SDSS J090334.92+502819.2 (hereafter SDSS J0903+5028), discovered in the Sloan 
Digital Sky Survey data (SDSS; (York et al. 2000)). There have been several gravitational lenses 
discovered previously in the SDSS data (Inada et al. 2003a; Inada et al. 2003b; Buries et al. 2003; 
Pindor et al. 2003b; Morgan et al. 2003), but this one is unusual in the way it was found. The stan- 
dard algorithm for selecting lens candidates in the SDSS involves looking for deviations from PSF 
profiles for spectroscopically confirmed quasars (Inada et al. 2003a; Pindor et al. 2003a). By con- 
trast, SDSS J0903+5028 was selected for follow-up based on the presence of z — 3.6 quasar features 
superimposed on the SDSS fiber spectrum of a z = 0.388 luminous red galaxy (LRG). Moreover, 
the SDSS image showed another source about 2 / .'5 from the spectroscopically targeted galaxy as 
well as the presence of surrounding galaxies with colors similar to the LRG. While this system 
may be unusual, it is not unprecedented: the well-known lens 2237+0305 (Huchra et al. 1985) was 
discovered serendipitously in a galaxy redshift survey, and such cases are expected (Kochanek 1992; 
Mortlock & Webster 2000b). 
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Follow-up spectroscopy at the ARC 3.5-meter telescope showed that both main image compo- 
nents contain flux from az~ 3.6 quasar with strong BAL-like associated absorption (Foltz et al. 1986), 
strongly suggesting that this is a lensed system. Subsequent r and i band imaging at the ARC 
3.5-meter telescope in better seeing revealed the lens geometry more clearly, as shown in Figure 1, 
and enabled us to model the lens system. Finally, higher signal-to-noise ratio, higher-resolution 
spectra were taken at the Keck II telescope. Taken together, these data make a firm case that 
SDSS J0903+5028 comprises two images of a high-redshift quasar lensed by a massive, red, fore- 
ground galaxy in a group. 

This paper is organized as follows. In Section 2 we describe the SDSS data on this object, 
discuss how it was selected for follow-up, and describe the spectroscopic and imaging data from the 
ARC 3.5m telescope and spectroscopic data from the Keck telescope. In Section 3, we fit a simple 
two-quasar + galaxy model to the ARC 3.5m images and extract positions and magnitudes for the 
three components. With this information, we fit a lens model, estimate the velocity dispersion of 
the lens galaxy, and study the quadrupole moment of the lensing potential. We also decompose 
the SDSS and Keck spectra into quasar and galaxy components and find flux ratios consistent with 
the imaging data. We conclude in Section 4. 

2. Observations 

The data on SDSS J0903+5028 consist of the following. The object was observed in routine 
SDSS imaging in January 2001. Based on its colors and brightness, it was targeted for SDSS 
spectroscopy as a luminous red galaxy (LRG) and spectroscopic observations were taken in March 
2001. Due to the presence of both quasar and galaxy features in the SDSS spectrum, the object 
was included in a list of promising lens candidates for follow-up observations. Spectroscopy of the 
two main components on the ARC 3.5m telescope in October 2002 revealed that they both contain 
very similar quasar features in addition to galaxy spectral features. This observation was followed 
by higher quality spectroscopy with the Keck II telescope and deeper, better seeing-quality imaging 
data with the ARC 3.5m telescope. In the following, we describe each of these data sets in detail. 

2.1. SDSS Data 

The Sloan Digital Sky Survey is a wide-field photometric and spectroscopic survey being carried 
out by the Astrophysical Research Consortium (ARC) at the Apache Point Observatory in New 
Mexico (York et al. 2000). The SDSS multi-CCD camera (Gunn et al. 1998) will produce images 
for ~ 5 x 10 8 objects in five optical bands u, g, r, i, z to a detection limit of approximately r = 22.2. 
The photometric pipeline software is described in Lupton et al. (2001). The photometric calibration 
is described in Fukugita et al. (1996) and Smith et al. (2002) and the astrometric calibration in Pier 
et al. (2003). Galaxies, quasars, stars, and other sources identified in SDSS imaging are targeted 



-4 - 



for SDSS spectroscopy based on several selection criteria (Stoughton et al. 2002); for a description 
of the selection algorithm for luminous red galaxies (LRGs), see Eisenstein et al. (2001). 

SDSS J0903+5028 was imaged on 2001 January 26 with SDSS identifiers: run 2074, cam- 
era column 2, field 113. This field is included in the recently released SDSS Data Release 1 20 
(Abazajian et al. 2003). SDSS J0903+5028 is the close pair of objects near the center of the r 
band image shown in Figure 1. The Western member of the pair is a blend of galaxy and quasar 
light and corresponds to the combined object (s) labeled G, B in Figure 1. In the photometric 
reduction used for spectroscopic targeting (rerun 0), this object was assigned SDSS identification 
number id 229 in this field. In the current 'best' reduction (rerun 21), using an improved version 
of the photometric pipeline, the object has id 186. The Eastern member of the pair — component A 
in Figure 1 — was also identified by the photometric pipeline (id 228 in rerun and id 185 in rerun 
21, with the same run, camcol, and field numbers as above). It was identified as a faint galaxy, not 
a point source, most likely because of its close proximity to the B/G component. 

In Table 1 we provide SDSS photometric (Hogg et al. 2001; Smith et al. 2002) and astrometric 
(Pier et al. 2003) parameters for both components from the 'best' rerun 21. Since the seeing 
measured in bands u,g,r,i,z was 1'.'83, l'/73, l'/50, 1'.'40, l'/50, respectively, comparable to the 
separation of the two images, the A and G/B components were not fully deblended from each other 
by the SDSS photometric pipeline. As a result, one should be cautious about interpreting the SDSS 
photometric parameters for this object. More accurate astrometric and photometric information 
for the different components, based on subsequent ARC 3.5m imaging in better conditions, is given 
in Table 2. 

The G/B component was targeted for SDSS spectroscopy as a 'cut II' luminous red galaxy 
(LRG) (Eisenstein et al. 2001; see Blanton, et al. 2003 for a description of the spectroscopic tiling 
algorithm) and was likely boosted above the flux limit of the LRG sample by the addition of the 
blended quasar light. Its spectrum was taken on 2001 March 24 with SDSS identifiers: Plate 552, 
Fiber 221, MJD 51992. This plate was observed for 8 x 15-minute exposures, yielding combined 
spectra of somewhat higher S/N than is typical for the survey. The 3" spectroscopic fiber was 
centered on 09 h 03 m 34 s .92 + 50°28'19?2 (J2000); in Figure 1, this corresponds approximately to 
centering on the G component. The SDSS spectrum, shown in Figure 2, clearly shows absorption 
features of an early-type galaxy at redshift z = 0.388 (e.g., the Ca H and K lines at 5463 and 5510 
A), along with strong quasar emission lines with a peak C IV redshift of z = 3.584. 

The surrounding field in Figure 1 shows several fainter galaxies with colors similar to those 
of the G/B component of SDSS J0903+5028, suggesting the presence of a small galaxy group as- 
sociated with the LRG. We have applied a group finding algorithm which looks for a red-sequence in 
color-magnitude space for over-dense regions, the maxBCG algorithm (Annis et al. 2003; Bahcall et al. 2003). 
The algorithm does not find a cluster or group centered at the lens galaxy since it is not a local 
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maximum of the galaxy density on cluster length scales. However it does find two small nearby 
clusters, both about 5'.7 away and both with photometric redshifts of z=0.44 ± 0.03. The lens 
galaxy sits right between these two clusters and they are aligned north and south of the lens. The 
properties of these photometric clusters are summarized in Table 3. The actual physics of these 
two clusters will require a detailed spectroscopic study: they may be part of one bigger cluster or 
set of merging clusters that span the entire region including the compact group of galaxies sur- 
rounding the lens galaxy shown in Figure 1, but the precision of photometric redshifts does not 
allow definitive answers to these kind of questions. 



2.2. Selection for Follow-up 

SDSS J0903+5028 was recognized as a possible gravitational lens during routine testing of the 
spectroscopic outputs of the SDSS. There are two independent SDSS software pipelines developed 
for classifying spectra and assigning redshifts: spectrold (briefly described by Stoughton et al. 2002; 
for more detail, see SubbaRao et al. 2003), which uses both cross-correlation via Fourier transforms 
with a family of templates and emission line identification, and specBS (Schlegel 2003), which uses 
X 2 template fits in wavelength space. Significant discrepancies in redshifts and/or classifications 
between the two were examined by eye 21 . One common type of discrepancy arises when light from 
superposed objects falls within the 3" spectroscopic fiber, and the two pipelines make different 
choices about which object's redshift to report. In DR1, there were half a dozen galaxy-quasar 
superpositions at very different redshifts identified in this way, including SDSS J0903+5028. None of 
the others appear to be lenses. In the case of SDSS J0903+5028, spectrold returned a classification 
of Galaxy with a redshift of z = 0.388 at 94% confidence, while specBS returned a classification of 
Quasar with a redshift (albeit incorrect) of z = 1.788. 

In addition to comparison of the two pipelines, the spectrold pipeline also flags spectra which 
cross-correlate with two templates at substantially different redshifts at high confidence level. Such 
was the case with SDSS J0903+5028: spectrold reported a significant (80% confidence) z = 3.6 
cross-correlation peak for this spectrum with a quasar template. The relative confidence levels of 
the galaxy and quasar peaks are in line with expectation, given that the galaxy flux through the 
3" fiber is about twice that of the quasar (see below). 

One of the unusual features of this lens is the fact that the target was identified as a luminous 
red galaxy, not as a quasar: the lensing galaxy is brighter than the lensed quasar images (see 
Table 2). This system would therefore not be included in many optical searches for lensed quasars, 
because the brighter component has galaxy rather than quasar colors and because it was identified 
as an extended rather than a point source. On the other hand, in surveys that extend to faint 



21 This totaled only 1.7% of all spectra in SDSS Data Release 1, of which roughly half were of too low S/N to yield 
a meaningful redshift, usually correctly classified as "unknown" by both pipelines. 
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magnitudes, it is not completely surprising to find such objects. For example, for UV-excess selected 
quasars at z < 2.5, Kochanek (1991) estimated that in a few percent of three-image lenses, the lens 
galaxy flux will exceed that of the combined quasar light for surveys to m = 21; presumably this 
percentage is higher for multi-band surveys that include quasars to higher redshifts. Alternatively, 
the typical r-band flux for a spectroscopically targeted z = 3.6 quasar in the SDSS is r ~ 19.4 (PSF 
mag), while the typical r-band flux from a targeted z = 0.38 LRG is r ~ 18.8 (model mag). While 
this comparison is obviously biased by our target selection criteria, it is nevertheless suggestive. 



2.3. ARC Spectrum 

We conducted follow-up spectroscopy of SDSS J0903+5028 and other interesting lens can- 
didates on 2002 October 9 with the Astrophysical Research Consortium (ARC) 3.5m telescope 
at Apache Point, New Mexico using the Double Imaging Spectrograph 22 . This instrument has a 
dichroic at 5550A; the red and blue spectra combined have a usable wavelength coverage of about 
3700A to lOOOOA. We took a 22 minute spectrum with the slit aligned along the direction con- 
necting the two primary image (A and B/G) components. Although the two spectra were partially 
blended and of relatively low signal-to-noise ratio, it was clear after reduction that both spectra 
contained flux from high redshift quasars at the same redshift of z ~ 3.6. Both spectra also showed 
absorption features from the galaxy. Because the DIS data established a strong case for lensing but 
were not definitive, we subsequently re-observed this lens candidate at the W. M. Keck Observatory. 



2.4. Keck Spectrum 

We obtained a high-dispersion spectrum of both image components of SDSS J0903+5028 using 
the echelle spectrograph and imager (ESI; Epps & Miller 1998) on the Keck II telescope on the 
night of 2002 December 5; see Figure 3. Three slit orientations were used; here we report only 
the pair of spectra taken with the slit perpendicular to the axis separating the image pair, since 
these observations yielded the cleanest reductions. The night was clear, with 0?8 seeing. A 900s 
high-resolution spectrum was taken for each member of the pair through a l" slit in the echellette 
mode of ESI. In this mode, the spectral range of 3900A to llOOOA is covered in 10 spectral orders 
with a nearly constant dispersion of 11.4 km s -1 pixel -1 . Wavelength calibrations were performed 
with observations of a CuAr lamp. The spectrophotometric standard BD+28 4211 was observed for 
flux calibration. The data were reduced using a tailored set of IRAF 23 and IDL routines developed 
specifically for ESI data. The smoothed spectra are shown in Figure 3. The signal-to-noise ratio is 



22 DIS II see http://www.apo.nmsu.edu/Instruments/DIS/ 

23 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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~ 5 pixeP 1 in the raw spectra at 1450A in the rest-frame and ~ 29 pixeP 1 in the smoothed spectra. 
As with the SDSS spectrum, the peak C IV redshift is z = 3.584. However, given that there is 
associated absorption long-ward of the emission peak, this redshift is likely to be an underestimate. 
Using templates from Richards et al. (2002), which allow for the possibility that C IV emission 
is blueshifted with respect to systemic and also allowing for reddening of the spectrum, we find 
a best-fit redshift of z = 3.605 for the less contaminated 'A' component, which would place the 
associated C IV absorption at roughly the systemic redshift (instead of being infalling). Figure 3 
shows that the two components have remarkably similar spectra and consistent redshifts. Scaling 
the fainter Eastern ('A') component by a factor 1.3 leads to a good match with the brighter Western 
('G/B') spectrum. Furthermore the spectrum is by no means a typical quasar spectrum since it 
has irregular BAL-like troughs. The fact that both components have these same rare BAL troughs 
makes the lensing case very solid; the fact that there is clearly a galaxy between them makes the 
case practically certain. 



2.5. ARC Imaging 

From the SDSS imaging data, it was apparent that SDSS J0903+5028 does not simply comprise 
two point sources: the G/B source is extended and was tentatively interpreted as a possible super- 
position of the lens galaxy with a quasar point source. To further test the lensing hypothesis and to 
determine source positions and magnitudes for lens modeling, we obtained follow-up imaging data 
on 2002 November 13 with the ARC 3.5m telescope using SPIcam. SPIcam is a backside-illuminated 
SITe 2048 x 2048 CCD camera with 2A^m pixels and a plate scale of 0'. / 14 pixel -1 , giving a field of 
view of 4'78. Because of the small pixels, this camera can take advantage of very good seeing. As 
it turned out, the seeing was l'.'l, a significant improvement over the SDSS 1"5. Also, the longer 
exposure, co-added SPIcam images are about 1.8 magnitudes deeper in r than the corresponding 
SDSS image. We obtained four dithers in each of the SDSS r and i bands for a total exposure time 
of 20 minutes in each band. These images were de-biased, overscan-corrected, and flat-fielded in 
the usual manner with IRAF. Figure 1 shows a 35" by 35" region around SDSS J0903+5028 from 
the co-added r band SPIcam image. The small group of galaxies is evident. Based on modeling (see 
below), the objects labeled A and B were identified as the quasar images, while the object labeled 
G is the galaxy image, blended with quasar image B. We used SExtractor (Bertin & Arnouts 1996) 
to find objects in the co-added r and i images and matched these to the SDSS imaging catalog to 
obtain photometric zero-points and an accurate astrometric solution. 
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3. Analysis 

3.1. Modeling the Image 

In order to fit a lens model to the data we proceed to determine the positions and relative 
fluxes of the quasar and galaxy images. While ideally one would like higher resolution images for 
this purpose, we can in fact determine the configuration of this system quite confidently with just 
arc-second imaging. The co-added ARC 3.5m images in both SDSS r and i filters are used to fit 
for an image model. 

The top left panels in Figures 4 and 5 show the r and i band co-added SPIcam images of 
the 8" by 5" area around the lens. The object on the left (East, component A) is unresolved and 
is one of the quasar images. The object on the right (West) is resolved and it is evident from 
visual inspection that this is in fact bimodal, with a point source, the quasar, to the lower right 
(southwest) of the blended object centroid. 

This hypothesis can be tested by fitting the image to a simple parametric model and looking at 
the residuals. The simplest model consists of a two-image lens with the galaxy in between the two 
quasar images. There are some conditions that must be met for this image to be consistent with 
gravitational lensing. The quasars should have identical shapes, consistent with the local point- 
spread function (PSF), while the galaxy may be more extended. The three objects should have 
positions in the two bands that are statistically consistent, and the quasars should have nearly 
identical flux ratios. There are other conditions that relate the flux ratios and the three image 
positions that arise from the gravitational lens model; we address those additional constraints in 
the Section 3.3. 

We fit the surface brightness of all three objects as two-dimensional t-distributions, also known 
as Moffat profiles. A normalized, 2D t-distribution is given by 

where 5 = (x— /Lt) T £ _1 (x— fi), the vector \x is the image centroid, and X is the 2x2 symmetric matrix 
of moments which determines the shape of the elliptical isophotes. The free parameter v determines 
the logarithmic slope of the asymptotic profile. In the limit v — > oo, the surface brightness becomes 
a Gaussian, <fi(x) — > (27r) _1 |S| -1 / 2 exp(— 5/2); even for finite v, <p is approximately Gaussian near 
the centroid. 

Our PSFs are well fit by v = 2, so we fix v to this value for the two point sources. The galaxy 
is also reasonably well fit by a v = 2 profile (but with different moments), but it is better fit by a 
v = 1 profile, so we use the latter. We further require that the quasars have the same moments. 
This leaves 15 free parameters: 3 pairs of centroid coordinates, 3 fluxes, 3 PSF moments, and 3 
galaxy moments. The fits are done independently in each band. The best fit values are presented 
in Table 2. 
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The best fit dereddened magnitudes for the lens galaxy are r = 19.59±0.06 and i = 18.86±0.07, 
where we have included all errors from shot noise, calibration, and model degeneracy. We can use 
the measured redshift of 0.388 to calculate absolute magnitudes in both bands and then use the 
L — a relation (Faber & Jackson 1976) to estimate the galaxy velocity dispersion in each band. 
We use i^T-corrections from Bruzual & Chariot (1993) and correct for luminosity evolution using 
Bernardi et al. (2003) to arrive at M r = -22.26 and Mj = -22.74. The galaxy is therefore very 
luminous ;about 3L* in both bands. Using the L — a relations and the L — a scatter from Bernardi 
et al. (2003), we estimate the velocity dispersion as a r = 206 ± 53 km/s and a% = 213 ± 54 km/s. 
Typical lens galaxy velocity dispersions are 200-300 km/s, so these values are not unusual; they are 
also consistent with the velocity dispersion inferred from the lens model below. 

The top two panels of Figures 4 and 5 show the SPIcam data along with the best fit model. 
The middle panels show the best fit model separated into the quasars and galaxy. The lower panels 
show the residual image (image— model) and a contour plot with the relative positions and fluxes 
of the three components. One can see that the r and i data give visually consistent results. The 
best fit models have reduced x 2 of 0.99 in r and 1.00 in i, indicating that the model is a good fit to 
both bands. The inferred quasar flux ratios (B/A) in r and i are 0.483 ± 0.012 and 0.461 ± 0.021, 
consistent at the 1.3 a level. The quasar separations are 2'.'83 ± Of! 02 in r and 2? 80 =t 0'.'03 in i, 
consistent at about the 1 a level. 

The image model also yields a measurement of the ellipticity of the galaxy light. The uncor- 
rected model galaxy ellipticity is e = (1 — r 2 )/(l + r 2 ) = 0.12 in the r image and 0.19 in i, with 
position angles of 12.0 and 18.2 deg (East of the North-South axis) in the two bands; here, r = b/a 
is the ratio of semi-minor to semi-major axis of the surface brightness distribution. However, these 
numbers do not take into account the extent and anisotropy of the image PSF. We correct the 
galaxy shape measurement by subtracting the second moments of the local PSF from the second 
moments of the G model image. Using these deconvolved moments, the estimated corrected galaxy 
ellipticity is e = 0.27 in r and 0.32 in i; the corresponding position angles are 24.8 and 30.8 deg. 
The estimated error on the inferred ellipticity is about 0.1. 

While the image modeling above does not rule out more complicated lens configurations, it does 
show that this image is consistent with the simplest configuration of a two-image lens. We also note 
that we have applied this image fitting procedure to the lower signal-to-noise ratio SDSS images 
and to subsequent CFHT images (taken in better seeing but with more complex PSF structure) 
with very similar results. 

3.2. Modeling the Spectrum 

As with the imaging, we have also attempted to model the various spectra of SDSS J0903+5028 
as a sum of quasar and galaxy components, using Principle Component Analysis (PCA). A large 
number of redshifted SDSS quasar spectra are used to construct eigenspectra ej(A), which form 
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an orthonormal basis in terms of which any other quasar spectrum can be expanded, /qso(A) = 
Cjej(A). Similarly, a set of galaxy eigenspectra are constructed from many SDSS galaxy spectra. 
Spectra can be usefully classified by their coefficients Cj, provided they can be accurately recon- 
structed when the series is truncated at relatively small N. Three eigenspectra span the range of 
most galaxy types; for quasars, more components are needed. Here, we used 10 quasar and 10 
galaxy eigenspectra, constructed from samples of several thousand SDSS quasar spectra and about 
100,000 SDSS galaxy spectra (Yip et al. 2003). (The same galaxy eigenspectra are used in the 
SDSS spectrold spectroscopic pipeline to classify galaxies.) 

To decompose a spectrum containing both quasar and galaxy components, for which the two 
redshifts are known, we simply assume it can be modeled as a weighted sum of the galaxy and 
quasar eigenspectra, where the coefficients are determined by minimizing the x 2 of the reconstructed 
spectrum fit to the true spectrum. An example of this 20-parameter fit (hereafter called Model 
I) is shown in the top panel of Figure 6, which shows the best fit to the SDSS spectrum of the 
G/B component; for this fit, the quasar flux is about 25% of the galaxy flux summed over this 
wavelength range. Unfortunately, given the nature of this procedure, it is difficult to assign an 
error to this value. 

In addition to the model above, we experimented with two other models with fewer parameters. 
In Model II, instead of using 10 galaxy eigenspectra, we fixed the galaxy spectrum to have the shape 
of the average Luminous Red Galaxy spectrum constructed by co-adding a large number of LRG 
spectra (Eisenstein et al. 2003). This makes use of the information that LRG spectra are quite 
homogeneous and that the G component has colors typical of an LRG. An example is shown in 
the lower panel of Figure 6, which shows a decomposition of the Keck Western (G/B) spectrum 
using this model. As with the SDSS spectrum, the reconstructed quasar spectrum is a reasonable 
first approximation to the observed spectrum; not surprisingly, this procedure does not capture 
the BAL-like features, since the parent sample of SDSS quasar spectra used to produce the quasar 
eigenspectra did not include quasars with BAL features. The ratio of quasar to galaxy flux for this 
model is 42%; for comparison, Model I for this spectrum yields a quasar /galaxy flux ratio of 48% 
and yields a galaxy spectrum with the general spectral shape of an LRG. The differences between 
this spectral decomposition and that for the SDSS spectrum are not particularly troubling: the 
Keck spectrum has higher signal-to-noise ratio, and it is based on a narrow slit with rather different 
aperture from the SDSS fiber spectrum. On the other hand, the quasar /galaxy flux ratio for the 
Keck West spectrum model is in good agreement with that inferred from the r and i band imaging 
given in Table 2. 

In Model III, on the assumption that the Keck East spectrum has little contamination by the 
lensing galaxy, we fit the Keck West spectrum to a sum of Keck East and the LRG template or to 
Keck East plus 10 galaxy eigenspectra. This model generally gave poor or unphysical fits, consistent 
with the fact that Figure 3 appears to indicate that the Western component is somewhat bluer 
than the Eastern component. The latter result is somewhat surprising: given the lens geometry 
shown in Figure 4 and the results in Table 2, one would naively expect the Eastern component to 
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be less contaminated by the red lensing galaxy than the Western component, assuming the East 
and West spectral components correspond approximately to the A and B/G image components of 
Figure 1. This apparent discrepancy may be due in part to the spectral extraction algorithm, errors 
in relative spectrophotometric calibration, placement of the slit, intrinsic reddening of the quasar 
spectra, differential reddening in the galaxy, or quasar spectral variability on a timescale shorter 
than the time delay between the images. 

Finally, we also attempted to measure the galaxy velocity dispersion from the quasar-subtracted 
galaxy spectrum, but it was too contaminated by residual quasar absorption features to obtain a 
reliable result. 

3.3. Modeling the Lens 

To extract physical properties of the lens galaxy and its environment and to further test the 
lens hypothesis, we proceed to make lens models using the astrometry and photometry from the 
model analysis of the images. The uncertainties on the relative positions are (/.'02 in r and (//03 
in i. For the fluxes, we broaden the error bars to 10% to account for variability, microlensing, etc. 
(see Dalai & Kochanek 2002). We use standard isothermal lens models, because they are consistent 
with the observed properties of other individual lenses, lens statistics, and the dynamics and X-ray 
properties of elliptical galaxies (Fabbiano 1989; Kochanek 1993; Kochanek 1996; Maoz et al. 1993; 
Rix et al. 1997; Treu et al. 2002; Rusin et al. 2003). For the modeling we use Vt m = 0.3 and £l\ 
= 0.7, although these only affect the reported velocity dispersions and time delays, and then 
only at the few percent level. We use standard non-linear least-squares lens modeling techniques, 
implemented in the lensmodel software by Keeton (2001b). The data provide eight constraints: two 
each for two quasar image positions, two for the galaxy position, and two fluxes. A minimal model 
has eight parameters: the galaxy position (2) and mass (1), the ellipticity of the mass distribution 
or alternatively shear and its orientation angle (2), and the quasar source position (2) and flux (1). 
Even minimal models therefore have A^j f = 0, and so we are always able to find models that fit 
the data perfectly. Hence to estimate the uncertainties on the model parameters, we repeatedly 
add random noise to the 8 data points and refit to obtain a distribution of fitted parameter values. 

The fact that the quasar images and the galaxy are not collinear indicates a non-negligible 
quadrupole moment in the lensing potential, which may represent ellipticity in the lens galaxy 
and/or tidal shear from mass in the environment of the galaxy (Keeton et al. 1997). The presence 
of mass ellipticity might be expected because the deconvolved galaxy light is elliptical; moreover, 
the mass could be more flattened than the light. The presence of shear seems likely because of the 
surrounding galaxy clusters. To consider both possibilities, we first examine two simple models: 
(1) a singular isothermal ellipsoid (SIE) model, where the quadrupole is due entirely to ellipticity, 
and (2) a singular isothermal sphere (SIS) plus shear model, where the quadrupole is due entirely 
to tidal shear. 
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Both SIE and SIS+shear models can fit the lens exactly, with the parameters given in Ta- 
ble 4. Both sets of models seem reasonable: SIE models require a mass ellipticity e = 0.5-0.6, 
slightly larger than the ellipticity of the light (~ 0.3), while SIE+shear models require a shear 
strength 7 = 0.15-0.18, typical of lenses in group or cluster environments (Keeton et al. 1997; 
Kundic et al. 1997; Kneib, Cohen, & Hjorth 2000). There are small differences between the r- 
band and i-band models due to differences in the deconvolved positions of quasar B and the galaxy 
in the r and z-band data, but the differences are only at the la level. 

The models yield an Einstein radius of I'M, corresponding to a velocity dispersion of 250 ± 4 
km s _1 for the lens galaxy. This number is consistent with the velocity dispersion estimates made 
with the L — a relations in Section 3.1; the estimate from the lens model may be higher due to the 
surrounding group slightly enhancing the image angular separation. The implied total magnification 
of the system is a moderate factor of 3-4. The models also predict that the time delay between the 
images should be in the range 57-72 h~ x days. Because the predicted delay depends on the relative 
amounts of ellipticity and shear (Witt et al. 2000), the usefulness of this lens for Hubble parameter 
analyses will depend on how well the ellipticity and shear can be determined independently. 

It is interesting to note that in both SIE and SIS+shear models the quadrupole moment of 
the lensing potential is oriented almost exactly north-south, while the corrected galaxy light is 
inclined at ~ 30° . A misalignment of more than ~ 10° usually indicates that the lensing potential 
has both ellipticity and shear with different orientations (Keeton et al. 1998; Kochanek 2002). It 
is pointless to fit models with unconstrained ellipticity and shear to SDSS J0903+5028, because 
such models are under-constrained. However, analyses of other lenses suggest that it is reasonable 
to constrain the shape of the model mass distribution using the observed shape of the light distri- 
bution. The orientation angles of the mass and light are strongly correlated and typically aligned 
to within ~ 10°, even if there is no clear relationship between the ellipticities of the mass and 
light (Keeton et al. 1998; Kochanek 2002). Figure 7 shows results for SIE+shear models where we 
either fix the shape of the model mass distribution to that of the light (panel a) or just require that 
the mass distribution match the light within assumed uncertainties of 10° in orientation and 0.1 in 
ellipticity (panel b). In both cases, the constraint on the mass orientation provides an important 
lower limit on the shear strength. In other words, under the reasonable assumption that the mass 
distribution is aligned with the light distribution, the misalignment between the galaxy and the 
quadrupole moment of the lensing potential directly implies the presence of shear from the lens 
environment (the group or nearby clusters). Adopting a constraint on the mass ellipticity would 
then yield an upper limit on the shear, but this result is less reliable because there is no strong 
evidence that the mass ellipticity should match that of the light (Keeton et al. 1998). We also 
note that if the two nearby clusters found by the maxBCG algorithm are indeed separate spherical 
clusters, they would produce a shear with the requisite north-south orientation; however, given 
their relatively large angular separation from the lens galaxy, one would expect them to produce a 
combined shear of only a few percent. 

Thus, the lens models suggest but do not conclusively reveal that the group or clusters around 



-13- 



the lens galaxy play an important role in the lensing potential. The best way to test this hypothesis 
would be to obtain spectroscopy for galaxies in the field of the lens, to confirm the cluster(s) and 
identify members, and to measure the centroid and velocity dispersion of the cluster (s). Those 
quantities could be used to estimate the shear from the environment, and compared with the 
predicted shear strength 7 ~ 0.1-0.2 and orientation 6 1 = 136°-174° to test and further constrain 
the lens models. 

4. Conclusions 

We have identified a lensed quasar candidate, SDSS J0903+5028, based on the superposition 
of a z = 3.605 quasar and a z = 0.388 luminous red galaxy in an SDSS spectrum. Follow-up 
observations with the ARC 3.5-m and the Keck II telescope have confirmed that this is a two- 
image gravitational lens system, with image angular separation of 2 / .'8. The lens model is consistent 
with a massive galaxy with a velocity dispersion of 250 km sec -1 . The lens geometry indicates 
a quadrupolar lensing potential which can be generated by an elliptical galaxy mass distribution 
and/or tidal shear from what appears to be a group of galaxies surrounding the lens. The misalign- 
ment between the quadrupole and the galaxy light suggests that there is indeed significant shear 
from the environment. 
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Fig. 1. — SPIcam r-band image of area around SDSS J0903+5028. North is up and East is to the 
left. The scale of the image is 35" across, the pixel scale is (X'14/pixel and the seeing is l'.'l. The 
objects labeled A and B are the quasar images; the galaxy is labeled G and is blended with quasar 
B. These other galaxies may be a small group or part of two nearby clusters. 
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Fig. 2.— SDSS spectrum of SDSS J0903+5028 (smoothed by 9 pixels). The error spectrum (also 
smoothed by 9 pixels) is given by the dashed line. Dotted lines mark the centers of Lya and CIV 
emission for z = 3.584. The flux units are 10 -17 ergs/s/cm 2 /A. 
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Fig. 3. — Keck spectra of SDSS J0903+5028. The dotted vertical lines show the location of Lya and 
CIV at the peak CIV redshift of z = 3.584. The best fit redshift is instead z = 3.605. Top: Western 
(brighter 'B/G') component. Second: Eastern (fainter 'A') component. Third: Fainter spectrum 
times 1.3 over-plotted on brighter component. Bottom: Ratio of scaled fainter component (A) to 
brighter (B/G) component. The larger contamination of the brighter component by the lensing 
galaxy makes the absorption line strengths appear different in the two spectra, when in fact they 
are very similar. 
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Fig. 4.— Model fits for the SPIcam r-band image of SDSS J0903+5028. Top left: Image. Top right: 
best fit model to two point sources and one extended source. Middle: Best fit model quasar and 
galaxy surface brightnesses. Lower left: Residuals between best fit model and the image. Lower 
right: Surface brightness contours for the 3 model components. 
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Model fits for the SPIcam i-band image of SDSS J0903+5028. For legend, see Figure 4. 
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Fig. 6. — Top: SDSS spectrum of SDSS J0903+5028, decomposed into galaxy and quasar compo- 
nents: Spectrum (black), Galaxy (blue), quasar (green), sum of Galaxy and quasar components 
(red). Bottom: Keck Western spectrum of SDSS J0903+5028, decomposed into LRG and quasar 
components. 
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Fig. 7. — Results from SIE+shear lens models. The solid curves show models in which the position 
angle (PA) of the mass is constrained to match that of the light, while the dotted curves show 
models for which both the position angle and ellipticity of the mass are constrained to match those 
of the light. There are separate curves for models of the r and i-band data, (a) The mass, PA, and 
ellipticity are fixed, (b) The mass, PA, and ellipticity are free parameters but constrained by the 
light, with assumed uncertainties of 10° and 0.1, respectively. Note that the Xmin = ® oru y because 
N doi = 0. 
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Table 1. SDSS Photometry 



Object 


Mag 


a 


r 


i 


z 


A 


Petro 
PSF 


21.20 ±0.18 
21.78 ± 0.08 


19.53 ±0.12 
20.22 ±0.06 


19.16 ±0.11 
19.58 ±0.04 


19.14 ±0.34 
19.27 ±0.09 


B/G 


Petro 
Model 


20.85 ± 0.39 
21.33 ±0.10 


19.28 ±0.14 
19.26 ± 0.02 


18.73 ±0.18 
18.50 ±0.02 


18.38 ±0.22 
18.17 ±0.04 



Note. — griz Petrosian (1976), PSF, and model magnitudes returned by 
SDSS photometric pipeline reduction rerun 21 for run 2074, camcol 2, field 
113, object id 185 (A) and 186 (B/G). PSF magnitudes are appropriate for 
point sources (component A, to the extent it is deblended), while model mag- 
nitudes (Stoughton et al. 2002) are used to define colors for LRG targeting 
(Eisenstein et al. 2001) (component B/G). The model magnitude errors only 
include residuals from the model profile fit and are therefore artificially low. 
These are asinh magnitudes (Lupton, Gunn, & Szalay 1999) and are corrected 
for Galactic reddening according to the dust map of Schlegel, Finkbeiner, & 
Davis (1998). u magnitudes are not reported since no significant flux was mea- 
sured in this band. More accurate photometry based on deeper ARC 3.5m 
imaging is presented in Table 2 below. 
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Table 2. ARC Photometry: Model Results 



Object 




(J2000) 


r 


i 


r — i 


QSO A 
QSO B 
Galaxy 


09 h 03 m 35 s 
09 h 03 m 34 s 
09 h 03 m 34 s 


.132 + 50°28'2tf'21 
.877 + 50° 28' 18'/ 75 
.925 + 50°28 / 19 / /53 


19.99 ±0.01 
20.78 ±0.03 
19.59 ± 0.02 


19.43 ±0.01 
20.27 ±0.05 
18.86 ±0.04 


0.56 ±0.02 
0.51 ±0.06 
0.73 ±0.05 



Note. — Dereddened magnitudes. The astrometry is well calibrated to the SDSS 
astrometry and so the dominant errors simply come from the fitting routine and 
this error is about / /07. The photometry is also calibrated to the SDSS and has a 
systematic error in the zero points estimated at 0.06 in both bands. The relative 
photometry could in principle be better but, due to the degeneracy between QSO B 
and the galaxy, the error on the relative magnitudes are at about the same level. We 
conclude that the color difference between the two quasars is consistent with zero. 



Table 3. MaxBCG Clusters 



(J2000) 


angle 


Distance 


z 


N 9 al 


09 h 03 m 43 s + 50°32'58" 


5^7 


1.36 


0.44 


12 


09 h 03 m 47 s + 50°24'54" 


5^8 


1.39 


0.44 


14 



Note. - The two clusters near the lens galaxy. The 
first column gives the cluster center J2000 coordinates as 
reported by the maxBCG algorithm. The second column 
is the separation in arc-minutes of the cluster center from 
the lens galaxy. The third is the separation in Mpc at the 
indicated redshift. The fourth column is the photometric 
redshift as reported by maxBCG. The estimated errors on 
the maxBCG redshift estimates are typically 0.02 for low 
redshift and about 0.05 for these higher redshift clusters. 
The fifth column is N ga i, a richness measure returned by 
maxBCG, the estimated number of L* and brighter galaxies 
in the cluster. 



-27- 



Table 4. Lens Model Results 



Type 


Band 


Re (") 


e or 7 


9 e or e 1 (°) 


Mtot 


At (hr 1 days) 


SIE 


r 


1.42 ±0.03 


0.47 ±0.04 


4.6 ± 2.8 


3.62 ±0.19 


67.4 ± 2.8 




i 


1.38 ±0.04 


0.57 ±0.06 


2.0 ±3.0 


3.19 ±0.23 


72.0 ±4.2 


SIS+shear 


r 


1.44 ±0.03 


0.15 ±0.01 


3.1 ±3.0 


4.13 ±0.21 


57.0 ±2.2 




i 


1.43 ± 0.04 


0.18 ±0.02 


1.0 ±3.2 


3.76 ± 0.26 


57.8 ±3.1 



Note. - - Col. 3 gives the Einstein radius. Col. 4-5 give the ellipticity e and position angle 
9 e for SIE models, or the shear 7 and position angle 9 1 for SIS±shear models. Col. 6 gives the 
total magnification. Col. 7 gives the predicted time delay (for a cosmology with £Im = 0.3 and 
^ A = 0.7). 



